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Abstract 

The quark-meson coupling (QMC) model is examined using Georgi's "naive 
dimensional analysis" . We argue that the QMC model is quite natural as an effective 
field theory for nuclei. 
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Recently we have developed an effective field theory for both infinite nuclear matter 
and finite nuclei, the quark-meson coupling (QMC) model. (The model originated with 
Guichon in 1988 [|TJ.) The QMC model may be viewed as an extension of QHD || || 
in which the nucleons still interact through the exchange of scalar and vector mesons. 
However, the mesons couple not to point-like nucleons but to confined quarks (in the 
bag). In studies of infinite nuclear matter Jl|, Q it was found that the extra degrees of 
freedom provided by the internal structure of the nucleon give an acceptable value for the 
incompressibility once the coupling constants are chosen to reproduce the correct satura- 
tion energy and density for symmetric nuclear matter. This is a significant improvement 
on QHD at the same level of sophistication. Furthermore, the model has been success- 
fully applied to finite nuclei within the Born-Oppenheimer approximation @, |6], [7|- It has 
been found that the QMC model can reproduce the properties of finite, closed-shell nuclei 
(from 12 C to 208 Pb) quite well. For this reason it also provides a very attractive framework 
within which to develop a microscopic understanding of the nuclear EMC effect [§, [J. 

In general, an effective field theory at low energy will contain an infinite number of 
interaction terms, which incorporate the compositeness of the low-energy degrees of free- 
dom, namely the hadrons |3], [TIJ, and it is then expected to involve numerous couplings 
which may be nonrenormalizable. Thus, one needs an organizing principle to make sen- 
sible calculations. 



Manohar and Georgi |Tl|] have proposed a systematic way to manage such compli- 
cated, effective field theories called "naive dimensional analysis" (NDA). NDA gives rules 
for assigning a coefficient of the appropriate size to any interaction term in an effec- 
tive lagrangian. (This NDA has been extended to an effective hadronic lagrangian for 
nuclei [fLOH .) After extracting the dimensional factors and some appropriate counting fac- 
tors using NDA, the remaining dimensionless coefficients are all assumed to be of order 
unity. This is the so-called naturalness assumption. If naturalness is valid, the effective 
lagrangian can be truncated at a given order with a reasonable bound on the truncation 
error for physical observables. Then we can control the effective lagrangian, at least at 
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the tree level. 

NDA has been already applied to QHD g [0J, the Nikolaus-Hoch-Madland (NHM) 
model 0, [13| and some nuclear phenomena |fT4| . (Concerning QHD, Furnstahl et al. [|H|] 



concluded that the relativistic Hartree approximation (RHA) leads to unnaturally large 
coefficients due to the treatment of the vacuum in terms of AA-pairs excitation. This 



means that the loop expansion in QHD does not work as well as one would like [JlOj, |15|.) 

Here we use NDA to see whether the QMC model gives natural coefficients. In brief, 
NDA tells us the following: for the strong interaction there are two relevant scales, namely, 
the pion-decay constant /„- (~ 93 MeV) and a larger scale, A ~ 1 GeV, which characterizes 
the mass scale of physics beyond the Goldstone bosons. The NDA rules indicate how those 
scales should appear in a given term in the effective lagrangian. The rules are: 

1. include a factor of 1//^ for each strongly interacting field, 

2. assign an overall normalization factor of (/^A) 2 , 

3. multiply by factors of 1/A to achieve dimension (mass) 4 , 

4. include appropriate counting factors, e.g. 1/n! for <p n (where is a meson field). 

Since the QMC lagrangian in the mean-field approximation (MFA) is given in terms of 
the nucleon (?/>), scalar (a) and vector (u and p) meson fields, we can scale a generic 
lagrangian component as 

— M^)\ff {ff ® « 

where V and r stand for a combination of Dirac matrices and isospin operators. The p 
meson field is denoted by b. The overall coupling constant ce mnp is dimensionless and of 
(9(1) if naturalness holds. 

As shown in Ref. ||, the basic result in the QMC model with MFA is that, in the 
scalar and vector meson fields, the nucleon behaves essentially as a point-like particle with 
an effective mass m%, which is given by a (relativistic) quark model for the nucleon (like 
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a bag) and depends on the nuclear density through only the a field, moving in a vector 
potential generated by the vector mesons. Because of their Lorentz-vector character, in 
infinite nuclear matter the vector mesons provide no effect on the nucleon structure except 
for an overall phase in the quark wave function, which gives a shift in the nucleon energy. 
In an earlier version of the QMC model [|l], f|, [5], || , we considered the effect of the nuclear 
medium on the structure of the nucleon alone and froze the quark degrees of freedom in 
the mesons. We call this version QMC-I. 

However, it is true that the mesons are also built of quarks and anti-quarks, and 
that they may change their character in matter [|7], To incorporate the effect of 

meson structure in the QMC model in MFA, we suppose that the vector mesons are 
again described by a relativistic quark model with common scalar and vector mean- 
fields, like the nucleon in QMC-I. In this case the effective vector-meson mass in matter, 
m*(v = uj,p), will also depend on the a mean-field. 

The a meson itself is, however, not so readily represented by a simple quark model 
(like a bag), because it couples strongly to the pseudoscalar (2n) channel and a direct 
treatment of chiral symmetry in medium is important [fTTf] . On the other hand, many 



approaches, including the Nambu-Jona-Lasinio model |T7[ , the Walecka model |18| and 
Brown-Rho scaling [^] suggest that the a- meson mass in medium, m*, should be less 
than the free one, m a . We therefore parametrized it using a quadratic function of the 
scalar field (by hand) 



1 - a a (g a a) + b a (g a a) , (2) 



with g a a in MeV, where g a is the cx-nucleon coupling constant (in free space). In Ref. 0, 
we chose three parameter sets: (a a ; b a ) = (3.0, 5.0 and 7.5 xlO -4 MeV -1 ; 10, 5 and 
10 xl0~ 7 MeV -2 ), called A, B and C, respectively. The values for the sets A, B and 
C were determined so as to reduce the mass of the a-meson by about 2%, 7% and 10% 
(respectively) at saturation density. This version, involving the structure effects of both 
the nucleon and the mesons, was called QMC-II. 
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The QMC lagrangian in MFA is given by 0, |7| 



£-QMC — £free + 4tn + ^QMC) 



(3) 



where £/ ree and C em stand for the free lagrangian for the nucleon and mesons and the 
electromagnetic interaction, respectively, while Dqmc involves (strong) interaction terms. 
For QMC-I and QMC-II, D^ MC is respectively given by % 



pint. 

'"QMC-I 



[m N - m* N (a)) - g^ uj - g p J 70& 



'0. 



(4) 



and 
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(5) 



where g u and g p are the u- and p-nucleon coupling constants, respectively. The effective 
masses in the medium, m* (j = N, a, uj, p), depend only on the scalar field, and mj is the 
free mass. If we introduce the field-dependent a-nucleon coupling constant by 



rn 



N 



m N - gAa)a, 



(6) 



the nucleon structure effect can be completely absorbed into the coupling constant g a (cr)- 
(Note that g a (0) = g a in (@).) 

To determine the cx-dependence of m*, we need a model for the structure of the hadron 
involved. For simplicity, we use the MIT bag model. (An alternative model can be found 
in Ref. ||.) Then the mass at nuclear density ps is found to take quite a simple form (for 
Pb ~ 3p 0) where p is the saturation density for normal nuclear matter) 0: 



rrij ~ mj 



no, , 
-j\.9*a) 



1 " %M 



(7) 



where Uq is the number of non-strange quarks in the hadron j and dj is a slope parameter 
for the hadron j, which is given by the second derivative of the mass with respect to the 
a field. 



There are some coupling constants to be determined: the a- and cu-nucleon coupling 
constants (g a and g u ) are fixed by fitting the binding energy (—15.7 MeV) at the correct 
saturation density (po = 0.15 fm~ 3 ) for symmetric nuclear matter. The p-nucleon coupling 
constant g p is used to reproduce the bulk symmetry energy, 35 MeV. The slope parameters 
and the coupling constants are listed in Table [I]. (In the actual calculation we take the bare 
quark mass and the bag radius of the free nucleon to be 5 MeV and 0.8 fm, respectively, 
and m n = 939 MeV, m a = 550 MeV, m u = 783 MeV and m p = 770 MeV. For details, 
see Ref. 0.) 

Table 1: Coupling constants (dimensionless) and slope parameters (in units of 10~ 4 
MeV" 1 ) for QMC-I and QMC-II (sets A - C). 



type 


9a 


9w 


9 P 


Oat 




a p 


QMC-I 


8.24 


8.17 


9.33 


8.80 






A 


6.95 


5.82 


8.35 


9.01 


8.63 


8.59 


B 


7.04 


6.31 


8.14 


8.98 


8.63 


8.58 


C 


6.94 


6.45 


8.07 


8.97 


8.63 


8.58 



Using (Q) and (|7|), the QMC-I lagrangian gives four interaction terms, while the QMC- 
II lagrangian offers 16 terms due to the internal structure of the nucleon and the mesons. 
Now we can calculate the dimensionless coefficients, ct mnp in (|l]), for each interaction term 
using NDA. Table |2] summarizes the interaction terms and the corresponding dimension- 
less coefficients. (We take A = m N in (p.) 

As seen in the table, the QMC-I model provides remarkably natural coupling constants, 
which lie in the range 0.5 - 1.0. In QMC-II, 14 or 15 of the 16 coupling constants can 
be regarded as natural. Only the large value of c 5oo for set C and c 6oo f° r sets A - C 
are unnatural. Since the coefficients, C0500 and C0600, are respectively proportional to a a b a 
and b\ (a CT and b a are defined by (^)), we can see that those unnaturally large numbers 
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Table 2: Interaction terms and corresponding (dimensionless) coupling constants. 



term 


Clinnp 


QMC-I 


A 


B 


C 




c 1100 


0.82 


0.69 


0.70 


0.69 


vpa 2 ^ 


c 1200 


-0.55 


-0.40 


-0.41 


-0.40 




ClOlO 


-0.81 


-0.58 


-0.63 


-0.64 


$ (l) 70^ 


ClOOl 


-0.92 


-0.83 


-0.81 


-0.80 


a 3 


C0300 





0.40 


0.67 


1.0 


a 4 


C0400 




-3.6 


-2.2 


-4.4 


a 5 


c 0500 




3.3 


2.9 


8.3 


a 6 


c 0600 




-22 


-5.7 


-22 


ouj 1 


C0120 




-0.77 


-0.78 


-0.76 


a 2 u 2 


C0220 




0.85 


0.87 


0.85 


a 3 u 2 


c 0320 




-0.71 


-0.73 


-0.70 


a A u 2 


C0420 




0.39 


0.41 


0.39 


ab 2 


C0102 




-0.75 


-0.76 


-0.75 


a 2 b 2 


C()202 




0.84 


0.86 


0.84 


a% 2 


C0302 




-0.70 


-0.73 


-0.70 


a% 2 


CQ402 




0.39 


0.41 


0.39 
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are due to the parametrization of the a mass in matter. In particular, a large value for 
the coefficient b a leads to unnaturally values for c 5oo and c 6oo- (This is the reason why 
C0500 and C0600 in set B are much closer to natural than in sets A and C. See also below 
(EDO We recall that the reduction of the a-meson mass in matter was the sole feature of 
the model which could not be calculated but was put in by hand. There is nothing within 
the QMC model which requires b a to be so large. Therefore, the QMC model itself can 
be regarded as a natural effective field theory for nuclei. 

Finally, we compare the QMC model with the NHM model JT2J and a new lagrangian, 
which was constructed recently by Furnstahl, Serot, Tang and Walecka [[3], |1(| (we call 
the latter the FSTW model). The NHM model was motivated by empirically based 
improvements to QHD , but using contact (zero-range) interactions to allow treatment 



of the Fock terms. It has 9 coupling constants, and 6 of them are natural [13]. On 
the other hand, the FSTW model was constructed in terms of nucleons, pions and the 
low-lying non-Goldstone bosons, and chiral symmetry is realized nonlinearly with a light 
scalar (a) meson included as a chiral singlet to describe the mid-range nucleon-nucleon 
attraction. This model has a total of 16 coupling constants, and they are almost all 
natural. In both cases, the coupling constants were determined so as to fit measured 
ground-state observables of several nuclei by a self-consistent procedure that solved the 
model equations for the nuclei simultaneously and minimized the difference between the 
measured and calculated quantities using a nonlinear least-squares adjustment algorithm. 
Therefore, the coupling constants were fixed entirely phenomenologically. 

This is quite different from the case of the QMC model. The latter has basically three 
coupling constants, which were determined to fit the saturation properties of nuclear 
matter (as mentioned after (0)), but the other coupling constants are automatically given 
through a model for the structure of the hadrons. Therefore, the physical meaning of 
the coupling constants is quite clear. In the QMC model, the meson mass (a, u, p etc) 
decreases in matter. However, in the FSTW model some parameter sets lead to an increase 
of the effective meson masses in matter, which seems unlikely from the point of view of 
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the recent discussions on this matter by many authors |J7|, [L8|, [19], |20], |21 



In conclusion, we have demonstrated that the QMC model is quite natural as an 
effective field theory for nuclei. Although we discussed the QMC model using a specific 
hadronic model, namely the MIT bag, the qualitative features we found here are expected 
to hold in any model in which the hadrons contain relativistic quarks. 



This work was supported by the Australian Research Council. 
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